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CLAIMS 

I claim: 

1. A method for downward extrapolation of pre-stack seismic data, comprising: 
selecting a set of prestack seismic data; 

determining a migration interval in the seismic data set; 
selecting a maximum error criterion for the migration interval; 

calculating a maximum relative error in phase as a function of frequency, propagation angle, 
and the relative variation in velocity in the migration interval; 
comparing the maximum relative error in phase to the maximum error criterion; and 
determining the type of extrapolation to use in the migration interval from the comparison of 
the maximum relative error in phase to the maximum error criterion. 

2. The method of claim 1, wherein the type of extrapolation is selected from a set comprising 
Gazdag phase-shift extrapolation, split-shift Fourier extrapolation, and implicit finite difference 
extrapolation. 

3. The method of claim 1, wherein the step of calculating a maximum relative error in phase 
comprises: 

selecting a frequency co in the migration interval; 

identifying a maximum wavenumber k max in the migration interval; 

identifying a minimum velocity c min in the migration interval; 

determining a maximum propagation angle 9 max in the migration interval, using the frequency 
co, the maximum wavenumber k max , and the minimum velocity c min \ 
selecting a reference velocity cq in the migration interval; 

calculating a maximum relative error in phase A& max , using the maximum propagation angle 
9 max and the reference velocity Co. 

4. The method of claim 3, wherein the maximum wavenumber k max is identified to include 95% 
of the power within a circle of radius k=k max . 
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5. The method of claim 3, wherein the minimum velocity c min is identified by scanning the 
velocity slice for the next migration step to identify a minimum velocity value. 

6. The method of claim 3, wherein the maximum propagation angle 0 max is determined by the 
equation: 
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where k max is the maximum wavenumber, c min is the minimum velocity, and co is the frequency. 



7. The method of claim 3, wherein the reference velocity cq is determined by the equation: 
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where c min is the minimum velocity and 0 max is the maximum propagation angle. 

8. The method of claim 3, wherein the maximum relative error in phase A0 max for split-step 
Fourier extrapolation is determined using the equation: 
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where & is phase, A& 'vs phase error over the migration interval, 0is the propagation angle with 
respect to the vertical, c(x, y) is the velocity, and c 0 is the reference velocity. 



27 



